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We study the evolution of perturbations during the domination and decay of a 
massive particle species whose mass and decay rate are allowed to depend on the 
expectation value of a light scalar field. We specialize in the case where the light 
field is slow-rolling, showing that during a phase of inhomogeneous mass-domination 
and decay the isocurvature perturbation of the light field is converted into a curva- 
ture perturbation with an efficiency which is nine times larger than when the mass 
is fixed. We derive a condition on the annihilation cross section and on the decay 
rate for the domination of the massive particles and we show that standard model 
particles cannot dominate the universe before nucleosynthesis. We also compare this 
mechanism with the curvaton model. Finally, observational signatures are discussed. 
A cold dark matter isocurvature mode can be generated if the dark matter is pro- 
duced out of equilibrium by both the inflaton and the massive particle species decay. 
Non-Gaussianities are present: they are chi-square deviations. However, they might 
be too small to be observable. 



I. INTRODUCTION 

Strongly supported by the recent observations of the Wilkinson Microwave Anisotropy 
Probe (WMAP) satellite 0,0,330, inflation has now become the dominant contender 
for generating adiabatic density perturbations with an almost flat spectrum. In the standard 
picture, the observed large scale density perturbations are due to the fluctuations of the 
inflaton field, which are created during a period of accelerated expansion by amplification 
of quantum vacuum fluctuations 01 . 

Recently, alternative mechanisms for generating density perturbations after inflation have 
been proposed. They all assume that the early universe is filled with at least one light scalar 
field whose energy density is negligible during inflation. Fluctuations in the field are ampli- 
fied during the inflationary phase with a quasi scale-invariant spectrum and their amplitude 
fixed by the energy scale of inflation, ^0 ~ where is the Hubble parameter at the 
horizon crossing. If the vacuum expectation value (VEV) of the light field (p during inflation 
is small, its density perturbation 5p^/ ~ 50/0 is larger than the perturbations generated 
during inflation. However, since the light field is sub-dominant, its perturbation is initially 
of the isocurvature type. Later, it can be converted into a curvature perturbation. For a 
successful conversion one of these two ingredients are necessary: (1) the field must come 
(close) to dominate the universe, or (2) it must induce fluctuations in a second component 
which eventually comes to dominate the universe. 

In case (1) the scalar field is called curvaton. In the curvaton scenario 0, 3, 
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the isocurvature perturbation stored in the curvaton field is transformed into a curvature 
perturbation during a phase in which oscillates at the bottom of its potential behaving as 
a dust-like component and thus dominating over the radiation. 

Case (2) has been inspired by the idea that coupling constants and masses of particles 
during the early universe may depend on the VEV of some light scalar field. This idea is 
motivated by supersymmetric theories and theories inspired by superstrings where coupling 
'constants' and masses of particles are usually functions of the scalar fields of the theory. An 
interesting proposal is that large scale perturbations could be generated from the fiuctuations 
of the infiaton coupling to ordinary matte r H^ . 13l |. These can be converted into curvature 



perturbations during the reheating phase fl4] when the infiaton decays into radiation and 
reheats different patches of the universe with different temperatures and energy densities 
(see however This mechanism of conversion has been called inhomogeneous reheating 

in An extension and generalization of this mechanism is that the masses of particles 
produced during the reheating are allowed to vary 16] . If sufficiently long-lived, the massive 
particles can eventually dominate the universe. Due to the fiuctuations in the masses the 
mass-domination process becomes inhomogeneous and can convert the density perturbation 
of (f) into a curvature perturbation. 

Some authors jl?, 3, 3 (^^^ ^^^^ 0]) have studied the evolution of large scale pertur- 
bations during the decay of the infiaton in the case of a fiuctuating decay rate. Two of these 
groups [l7i lisl make extensive use of the formalism introduced in Ref. 21 1 (see however 



also 122, |23|) for interacting fiuids. Indeed, on several occasions the light scalar field can 



be treated as a fiuid and the formalism of Ref. [21J consistently applied. In this paper we 
extend the analysis of jl?! . 3] and we include the case where the inhomogeneous reheating 
is due to the decay of massive particles whose mass and d ecay rate depend on the VEV of 
a light field. In the limit of a constant mass the results of [l7|, |l8| are recovered. 

We also derive the condition for the domination of the massive particles. Since at early 
times they are in thermal equilibrium, their abundance depends on their annihilation cross 
section. As we shall see the condition for the domination does not depend on their mass 
but only on their annihilation cross section and decay rate. Our results are then compared 
with those obtained for the curvaton model and some of the observational consequences are 
discussed. If the infiaton and/or the massive particles decay out of equilibrium isocurvature 
perturbations can be generated. The inhomogeneous infiaton decay and mass-domination 
can also lead to non-Gaussianities in the adiabatic spectrum of perturbations. These can 
lead to observational signatures in the cosmic microwave background (CMB) that can be 
possibly observed in future experiments. 

The paper is organized as follows. In the next section we model the inhomogeneous mass- 
domination mechanism by writing down the conservation laws for the radiation, the massive 
particles and the light field, and we derive the coupling between these species. In Section 
III we write the background and perturbation equations obtained from the conservation 
laws. These equations are completely general. They can be used for the inhomogeneous 
infiaton decay as for the inhomogeneous mass-domination model. In Section IV we con- 
centrate on the limit that the light field is in slow-roll and we discuss the efficiency of the 
inhomogeneous infiaton decay and mass-domination models. We also compare them to the 
curvaton model. Finally, in Section V we discuss the observable signatures of the inho- 
mogeneous mass-domination i.e., the possible presence of isocurvature perturbations and 
non-Gaussianities. In the final section we draw our conclusions. 
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II. MODELING THE COUPLINGS 

As mentioned in the Introduction, we want to derive a set of equations that, depending 
on the initial conditions, can describe two different physical situations: (A) inhomogeneous 
inflaton decay or (B) inhomogeneous mass-domination and decay. In order to describe these 
two cases we consider a universe containing three components: a radiation fluid 7, a non- 
relativistic fluid ip, and a light scalar field (p. According to the different physical situations 
we have: 

(A) Inhomogeneous inflaton decay: At the end of inflation the universe is dominated by 
the inflaton. This coherently oscillates at the bottom of its potential behaving as a 
non-relativistic fluid ip. The decay rate of the inflaton, F, depends on 0. 

(B) Inhomogeneous mass- domination: After inflation and reheating, the universe is filled 
with radiation 7 and with a species ip of non-relativistic massive particles, whose mass 
M and decay rate T depend on (p. 

As we shall see, the equations derived for case (B) can be applied to case (A) in the limit that 
M does not depend on (p. Therefore, in the following we concentrate on the inhomogeneous 
mass-domination and decay. 

The universe is initially dominated by a radiation bath 7 with temperature T and by 
a massive particle species with mass M > T, such that it can be described as a non- 
relativistic fluid ip. If the massive particles froze out at a temperature such that M/T was 
not much larger than 1, then the species can have a significant relic abundance and, if 
sufficiently long-lived, can eventually dominate the universe before decaying. The condition 
for the domination can be obtained by requiring the decay rate to be smaller than the 
Hubble parameter at the moment when ip decays. The latter corresponds approximately to 
~ M^/Mpi if we assume that ip starts to dominate as soon as it becomes non-relativistic, 
T ~ M. Thus we find Q, 

r < MVMpi, (1) 

where Mpi = (SttG*)"^/^ is the reduced Planck mass. A more precise condition for the domi- 
nation is derived in Sec. IV. B and depends only on the annihilation cross section and decay 
rate of the massive particles. When the ■j/'-particles decay into 7 they release a considerable 
amount of entropy and reheat the universe once again. During the inhomogeneous mass- 
domination and reheating, fluctuations in M = M{(p) and F = T{(p) can be converted into 
curvature perturbations. 

Variations of masses and coupling constants can be easily motivated in the context of 
string theory where compactified extra dimensions emerge as massless scalar fields, called 
moduli, in the effective four dimensional theory. These moduli generally couple directly to 
matter leading to variations of masses and fundamental constants 23|. That the mass of 



particles may depend on a light field has been discussed in several places in the literature, 
especially in the context of interacting dark matter and dark energy (see e.g., pil for a 
review) and more recently in models with variation of the fine structure constant 26] and 
in the chameleon cosmology j23|. A field dependent mass leads to a non-trivial coupling 
between the field and the non-relativistic fluid made of massive particles. The form of this 
coupling can be derived by considering the non-relativistic fluid as a classical gas of point- 
like particles with mass M{(p) and action S = — J Mds |23|. In this section we derive the 
conservation equations of our three fluids system using the conservation of the total energy. 
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Let Tl2\ T'^\ and be the energy-momentum tensors of the three components 7, ij)^ 
and 0. In particular we write the energy-momentum tensor of the non-relativistic massive 
species as 

T^^^ = Pi^u^Uy = Mn^u^Uy, (2) 

where is the energy density and the number density of ip. From general covariance 
we require that the sum of the three energy-momentum tensors is conserved, 

= V^(Ty)'^ + Ti^)'^ + TW'') (3) 
= V^Ti^)^' + V^Tl^^^ + ^T^^ + MV,{n^u^)uy + Mn^u^V ,Uy. (4) 

The last term of the right hand side of this equation vanishes only if the massive particles 
follow a geodesic, which is not the case if they interact with other fluid components. We 
project Eq. (j^ on the unit vector u'^ and we define the decay rate of ip as the rate of change 
of the number of ^/'-particles, 

V^{n^u^) = -Tn^. (5) 

By assuming that ip decays only into 7, from Eq. we obtain three coupled conservation 
equations, 

M-V^Ti^)^ = -{CmO.cPu'^ - r)p^, (6) 
^.y^TW^ = -Tp^, (7) 
^ Cm^^^m^, (8) 

where we have defined the mass coupling function Cm, 

CM ^ (9) 

For later convenience, we also define the decay rate coupling function Cr, 

CK.).^^. (10, 

Equations (jUHHI) will be used to derive the background and perturbation evolution equations 
describing the inhomogeneous mass-domination. The perturbation equations for the coupled 
field-fluid system (F = 0) in the case of a linear coupling (Cm = constant) have been also 
discussed in ji^ in the context of coupled quintessence. Finally, if Cm = 0, Eqs. (jHHHl) 
describe the inhomogeneous infiaton decay, as discussed in 



t maiiy, 

Ii3l3. 



III. COSMOLOGICAL PERTURBATIONS 



Here we write down and discuss the background and perturbation equations derived from 
Eqs. (HHD- 
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A. Background equations 

Consider a Friedmann-Lemaitre universe, with metric ds^ = —dt^ + a?{t)dx^, governed 
by the Friedmann equations, 

where p and P are the total energy density and pressure of the universe, and are the 
energy density of the 7 and species, respectively, and H is the Hubble parameter, defined 
as if = d/a. 

The evolution equation governing the total energy density and ensuring the conservation 
of the total energy is 

p = ~H (3p^ + 4p^ + 302) . (13) 

In particular we have three coupled conservation equations which can be derived from 
Eqs. (PIHl)- They read 

p^ = -Cm^ + T)p^, (14) 

p^ = -AHp^f + Fp^, (15) 

P0 = -?>H(t? - CM4>Pi'- (16) 

In order to solve for we need the field evolution equation, 

4> + 3H<P + V^ = -CmP^, (17) 

where V(j, is the derivative of the potential of 0, \^ = V{(j)), with respect to the field 0. 
As first introduced in ji^] we define the energy transfer functions Qa (with a = ip,'j,(j)) 

as 

Qa = Pa + 3H{p^ + P^). (18) 

From Eqs. ()14lll6p we have 

Qi. = (CM0-r)p^, (19) 

= Fp^, (20) 
= -CM(ppi>- (21) 



B. Perturbation equations 

Here we discuss the perturbation equations derived from Eqs. (jHHHl)- We describe scalar 
perturbations in the metric with the line element 

ds"^ = -(1 + 2^)df + a^{l - 2'^)6ijdx'dx\ (22) 

where $ and correspond to the Bardeen potentials in longitudinal gauge. In the absence 
of anisotropic stress perturbation $ = \E'. In order to perturb the energy- momentum ten- 
sor of the three components -0, 7, and 0, we introduce the energy density and pressure 
perturbations, 6pa and 6 Pa- 



6 



For each fluid component a one can introduce the gauge invariant curvature perturbation 
on the uniform a-energy density hypersurface defined as 

= _^ _ if^. (23) 

pa 

The Ca's remain constant on large scales only for adiabatic perturbations and in any fiuid 
whose energy- momentum tensor is locally conserved: n''Tjf'>^ = [s^- In our case Qa 7^ 
and this is not the case. The total uniform curvature perturbation, introduced by Bardeen 
|29| and Bardeen, Steinhardt and Turner [s^, is given by a weighted sum of the individual 
uniform curvature perturbations, 

C = - = E -Cc-, (24) 

P a P 

and the relative entropy perturbation between two components a and is given by the 



difference between the two uniform curvature perturbations (a and [21 



S^p = 3iCa-Cp). (25) 

According to Ref . Ill] , at large scales (neglecting spatial gradients) the time evolution of 
the curvature perturbation can be written as 

C = ^SP^.d, (26) 

p + P 

where 5Pnad is the total non-adiabatic pressure perturbation, given by the sum of the intrin- 
sic non-adiabatic pressure perturbation of each component and the relative non-adiabatic 
pressure perturbation, 

SPn.d = T.^Pintr,a + SPreU (27) 

a 

where the intrinsic non-adiabatic pressure perturbation of the species a is given as 

SPintr,a = 5P^ - ^Jp^, (28) 

and the relative non-adiabatic pressure perturbation depends on the uniform curvature per- 
turbations, 

5Pre^ = E ^(C - C.) = -T^.HMcl - 4)S^,. (29) 

Here = Pa/ Pa is the adiabatic speed of sound of the species a. At large scales, the 
evolution equation for each individual uniform curvature perturbation C,a is given by [2l| 

Pa 

The Ca's are sourced by three gauge invariant terms, the intrinsic and the relative non- 
adiabatic energy transfer functions, and the intrinsic non-adiabatic pressure perturbation. 

We now describe, one by one, the three terms on the right hand side of Eq. (|3U|). The non- 
adiabatic pressure perturbation (5Pintr a is defined in Eq. . For a fiuid whose parameter 
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of state Wa = Pa/ Pa IS constaiit, = Wa and (5Pintr,Q = 0. This is the case for the non- 
relativistic species and for the radiation 7, 

'^-Pintr.V' = '^-Pintr,7 = 0. (31) 

However, a scalar field generically has a non- vanishing non-adiabatic pressure perturbation 

13, 



5Pi„tr,, = 2Vj -2V^ (q, + IC^I , (32) 



where we have defined the gauge invariant variable 



Q^^6<P+j^^, (33) 



sometimes referred to as the Sasaki and Mukhanov variable j33|. If the field is the only - or 
the dominant - component of the universe, the intrinsic non-adiabatic pressure perturbation 
is negligible on large scales (see e.g., j3). However, in our case the contributes only to a 
small amount of the total energy and SPintr,<t, cannot in general be n egle cted. 
The intrinsic non-adiabatic energy transfer function is defined as j2l| 

SQmtv,a = 5Qa - ^6pa, (34) 
Pa 

where 5Qa is the perturbation of the energy transfer function of Eqs. (jl9H2H) . It auto- 
matically vanishes if the energy transfer function is only a function of the local energy 
density pa - when neither the decay rate nor the mass of the 'j/'-particles depend on the light 
field (f). For our three components ip, 7, and 0, we have 

^Qintr,^ = P^ I (cMi + Cm - ^) + jjQ^ + Cm (q<^ + 1:^) I , (35) 

^Qintr,7 = -P^ |- J (q<P + ^(7^ + ^"^'^'^T'l ' (36) 

SQmtT,(j, = —pip I Cm 7^77— + \Cm—^ — Cm] \ Q(t> + '7tC<H \ ■ (37) 
[ 6H \ (j) J \ ^ J ) 

The relative non-adiabatic energy transfer function is due to the presence of relative entropy 
perturbation and is given by 

^Q^el,a = ^ (C - Ca) = PP^^P- (38) 

It vanishes if the background transfer function = 0. 

Now we have all the ingredients to write the evolution equation of the total curvature 
perturbation ( on large scales, which is obtained from Eq. (jSHl), 



'^fe*^4cl (39) 
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and the three large scale evolution equations of the uniform curvature perturbations (a, 
which are derived from Eq. (jSDI), 







H 



Pj [ 3H 
Pj^ \ Cm4> 




(40) 
(41) 



( p± 






2p) 



Zp 



<t>i 





(42) 



where we have used F = CpF^. In order to close the system, we need the evolution equation 
for the Sasaki- Mukhanov variable on large scales, namely. 



H 



0C + -Q<i- 



Q6 + 




(43) 



Equations ()H9ll4Hj) are five coupled first order differential equations - one of them is redun- 
dant - which can be solved in order to study the evolution of perturbations during the 
inhomogeneous reheating. We did not make any assumption for the intrinsic non-adiabatic 
pressure perturbation of the field ^Pintr,0, which can be found from and (^p [see Eq. (IH^]. 
as well as for the coupling functions Cm and Cr- Thus these equations hold for any type of 
functional dependency of the decay rate and mass on the light field. 

Some comments are in order here. The evolution equations of the relative curvature 
perturbations and are sourced by terms proportional to p^, although has an extra 
term proportional to its intrinsic non-adiabatic pressure. After the decay of tp, is constant 
whereas ( and evolve only due to the intrinsic entropy of the light field. If the coupling 
function Cm vanishes, Eqs. ()H9II4H|1 describe the inhomogeneous inflaton decay, where ijj 
represents the oscillating inflaton and 7 the radiation produced by the decay and reheating. 



IV. SLOW-ROLL LIMIT 

As an application, we consider the limit where the scalar field is slow-rolling, 

~ 0, 4>/{H^) ~ 0, (44) 

and we assume that these conditions are maintained during all the period of ip domination 
and decay. 

In the slow-roll limit the scalar field varies very slowly in a flat potential and its energy 
density is always negligible. In this limit ~ — 1 and its intrinsic non-adiabatic pressure 
perturbation vanishes from Eq. 
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Although the variable diverges in this limit, the Sasaki-Mukhanov variable can still be 
used to discuss the scalar field perturbation. Equations ()H9ll4Hj) take a very simple form, 



C = J {pAC^ - C) + ^CmHp^Q^} , (46) 

= (^(C, - C^) + CrHaA , (47) 
Pi> ["^P ) 

= 0. (49) 

Also the background equations take a very simple form. In order to solve them numeri- 
cally, it is convenient to work in terms of dimensionless quantities, the density parameters 
Qa = Pa/ P, with + = 1 (fl^ ~ 0), and the dimensionless reduced decay rate j2ll |. 

j-r^. (50) 

which varies monotonically from to 1. In the slow-roll limit of Eq. ()44j) the background 
equations ()14II16|) can then be written as an autonomous system of first order differential 
equations, 

% = (^7 - (51) 
% = (^Y^ - (52) 

g' = lg{l-g){A~Q^), (53) 

where the prime denotes differentiation with respect to the number of e-foldings iV = Ina. 
Since these equations are subject to the constraint fl^ + fl^ = 1 there are only two indepen- 
dent dynamical equations whose solutions follow trajectories in the compact two-dimensional 
phase-plane {g,Q^). One can find a detailed analysis of this system, applied to the study 
of the curvaton model, in Ref. [2l|, where it is shown that close to the origin fi^ ~ Ri^^g^/'^, 
with 

i?in = , (54) 

in 

where the initial conditions for this system are set at t = tin such that g '^1. The initial value 
i?in determines which trajectory is followed in the two-dimensional phase-plane. For /2in > 1, 
the massive species ip comes to dominate the universe before decaying - compare with Fig. |21 
Indeed, the physical interpretation of -Rin is straightforward: If initially F ^ HQ'^, the decay 
is almost instantaneous and ip does not have the time to dominate. 

The perturbation equations, Eqs. (j46H49|) . can be written in terms of the dimensionless 
background quantities defined above and acquire a simple form 

C'-4^{(^)k*-C)-^C:„qJ. (55) 
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= 3^ {l^^ " ^^^^^ ' ~ ^"^4 ' ^^^^ 
Q; = 0, (57) 

where we have ehminated the variable which is redundant. When Cm = we recover 
the equations studied in with %jj representing the inflaton during its coherent oscillation. 
When also Cp = we recover the equations studied in with ip representing the curvaton. 

To calculate the final curvature perturbation on large scales, we start with initial condi- 
tions ai g <^ 1. By using the slow-roll condition (01)), from Eqs. and ()14H16p the initial 
total curvature perturbation is given by 



Cin — p^(3^^vC»/' + ^^-yC-r ~ CM^tpQi 



(58) 



and in the two physical situations that we discuss in this section it is negligible. Hence we 
shall consider vanishing initial perturbations, Cm — C7,in — Cv'.in — 0. Then we numerically 
solve the system of Eqs. ()55ll57p and we evaluate the perturbation variables for g ^ 1. The 
late time solutions approach a fixed point attractor, 

C = C7 = -r{arCr + aAfC'A/)-% = -ra„-% (59) 
= -r ar— + , (60) 



r M 

= -r(ar + l/2)-^-raA/-^, (61) 

where r < 1 is a function of i?in, whereas ar and gm are constant which can be determined 
numerically. The term ap -l- 1/2 in Eq. ()6H1 comes from the fixed point ('^ = for = 
0, (yf = 1 in Eq. (jSUj) . An example of the evolution of perturbations as a function of g is given 
in Fig. HJ 

The behavior of r as a function of i?in is illustrated in Fig. |21and it is obtained by solving 
numerically Eqs. ()55II57|) . For large i?in, r = 1 and the efficiency of the mechanism of 
conversion of the density perturbation of into a curvature perturbation is simply given by 
a^. For small Rm, r 1 and the efficiency is considerably reduced. In Fig. |21the maximum 
value reached by before the decay is also shown. For large R^^ one can approximate r by 
this value. Note that r as defined here is the same as the one defined in |lO| and computed 



m 
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for the curvaton model. We numerically checked that this is the case j46|. 
We can represent the integration of Eqs. ()55II57|1 as a transfer matrix acting on the initial 
perturbations, 

(62) 

Since the energy density of the light field is always negligible, its perturbation represents 
an entropy perturbation which is converted into a curvature perturbation if and r are 
non-zero. Our task is now to estimate the efficiency parameters ar and aM- Below we 
discuss two different physical situations. 
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FIG. 1: Evolution of C, Ci/)) ^'^d (j normalized to the light field perturbation —Q^/4>, as a function 
of g, for i?in = 10, where we have taken the simplest case where the dependence of V and M is 
linear in the field, M (x T (x (p. Although it is not visible on the fl')-plane, for g ^ 1 the 
uniform-density perturbations approach a fixed point attractor. 



A. Inhomogeneous inflaton decay 

Here we discuss case (A) as mentioned in Sec. II. The oscillating inflaton is described 
by the fluid ip and Cm = is used. Although fluctuating, the decay rate of the inflaton is 
time independent in the limit ()44|) . Initially, the inflaton dominates the universe, VL^ — > 1 
for (yf — i> 0, whereas the radiation is negligible; thus we take R\n 3> 1 and r = 1. The 
initial condition of C?/; is given by the initial perturbation of the inflaton field corresponding 
to vacuum fluctuations. According to [12] we assume that the density perturbation of the 
inflaton field is negligible. Therefore we have ^ \Cip,in\ = |Cin|, where we have used 

Eq. with Cm = and fi^ = for the last equality, and 0* is the VEV of at horizon 
crossing. Solving Eqs. ()55l | 57l) with R^^ ^ 1 (inflaton domination) we find ar = 1/6. This 
leads to the result found in [12e XL JA, 2M ; 

C = -g^, (63) 

valid on the spatially flat slices \1/ = 0. 



B. Inhomogeneous mass-domination and decay 

Here we discuss case (B) as mentioned in Sec. II. Now is the massive particle species 
and both the mass M and decay rate T depend on (p. The radiation 7 is the product of a 
previous reheating and it initially dominates the universe, fi^,in — 1. 

We start by discussing the condition for the massive particles to dominate the universe 
before decaying. Assuming that the massive species ip is sub-dominant when the initial 
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FIG. 2: The function r, parametrizing the conversion from an isocurvatm'e to a curvature pertur- 
bation, as a function of the initial parameter R[n (soHd line). This is compared to ri^^max (dashed 
line), the maximum value reached by fi^ before decaying as a function of i?in = 0^(i?/r)^/^|in. 



conditions are setup, the initial parameter can be written as 

T2 \i/2 



-Rir 



9^ 



1/4^ 



(64) 



where Qr is the number of relativistic degrees of freedom and n^/s is the reUc abundance of 
the ip particles when they freeze-out (s is the entropy density). In order to derive Eq. 
we have used 

s = ^grT\ (65) 



30^ ' 



45 



After freeze-out and when T <^ H,n^/ sis constant if g,- is constant (which we shall consider 
throughout). The massive species has time to dominate the universe if i?in > 1, which 
translates into 

M2 . /«,x2 



F < 



M, 



-9^ 



1/2 



PI 



(66) 



If n^/ s is order unity Eq. ((T)) is recovered. However, it is interesting to try to plug some 
numbers for the relic abundance in Eq. (j66|) . We make use of the^analytic approximation 
for the relic abundance of long-lived massive particles derived in |36] . At high temperature 
(T ^ M) oc T^, whereas at low temperature (T -C M) the density is Boltzmann 
suppressed, oc T^/^ exp(— M/T) so that if the particles freeze-out when T > M then 
the abundance becomes very small. The initial equilibrium abundance is maintained 
by annihilation of particles and antiparticles with cross section which we take to be 
independent of the energy of the particles. In this case the abundance at freeze-out is 
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r[GeV] 



FIG. 3: Values of F and {cav) for which R\a = 1. For values below this line the i/'-particles come 
to dominate the universe. We have taken — 100. 



where {ctav) is the thermal average of the total cross section times the relative velocity v. 
On using this relation, Eq. (jMj) becomes independent of the mass M, 

100 /Mpi\i/2 Mpi=^ 



This relation holds if the ■^/'-particles are subdominant at the freeze-out. If a more detailed 
calculation is p erformed one can see that the mass dependence enters only via a logarithmic 
correction [36|. In Fig. El we show the values of F and {(Tav) for which R^^ = 1 and the 
^/'-particles come to dominate. On requiring that the massive particles decay before nucle- 
osynthesis i.e., F > (1 MeV)VMpi, we find that only for (aAv) < (100/^^)i/2(^10^GeV)-2 
does ip dominate the universe. In this case, the initial thermal equilibrium by annihilation of 
particles and antiparticles must be maintained by some gauge interaction much weaker than 
those of the standard model. This excludes ip as being made by standard model particles. 

Let us now discuss the perturbations. The initial perturbation of the radiation is left 
over by infiation and it is negligible. The species ip is initially in thermal equilibrium with 
7, CiIjM = C7,in, and Eq. dHH)) implies Q^/cj)*. > C7,m = Cm- For the value of um we find 
numerically aM = 4/3. Thus, in general for a 0-dependent mass and decay rate 

C = -rfi^ + l^V (69) 
^ V6 F 3 M y ^ ' 

If the massive particles dominate the universe before decaying, i.e., if i?in ^ 1, r = 1. For 
M = constant we recover the result of the inhomogeneous infiaton decay discussed in the 
preceding section. If the mass depends on while F does not, we find 

C = (70) 
^ 3 M ^ ^ 



Finally, if M oc F we find 
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a result nine times larger than the one obtained in the case where only T is fluctuating. 

In Ref. fld\ the perturbation generated by a varying mass is derived with an analytic 
argument which yields ( = {1/3){Qm/M), a result different from the one of Eq. (fTHjl . This 
difference is due to the fact that in our model, specifled by the conservation equations (0- 
IH}, the entropy perturbation stored in the light fleld is important. Indeed, when Cm 0, 
p<^ 7^ and the light field perturbation contributes to the relative non-adiabatic pressure 
6Prc\ [see Eq. (j^ ] and sources the evolution of the total curvature perturbation, as shown 
by Eq. 



C. Comparison with the curvaton model 

It is worth discussing here the curvaton model, a mechanism of generation of perturba- 
tions which has very similar properties as the model discussed here. The aim is to stress 
their similarities and compare their efficiency. The curvaton a is a scalar field which is prac- 
tically free during infiation and starts oscillating after infiation (but before nucleosynthesis) 
during the radiation era when T < H < rrifj, behaving as a non-relativistic fiuid. If during 
the oscillations it comes (close) to dominate the universe before decaying, its perturbation 
(a is converted directly into curvature perturbations, 

C = rC.. (72) 

We can write the parameter i?in in terms of relevant quantities. If we choose the initial time 
when the curvaton starts to oscillate, at if ~ nia, on using p^- ~ m^a^ we have, from the 
definition (j3D), 

If the decay rate is sufficiently smaller than the curvaton mass, the curvaton has the time to 
dominate the universe before decaying. The conditions for the curvaton domination during 
its oscillations are 

r 

— Mpi < a < Mpi, (74) 

where the last inequality ensures that the curvaton does not dominate before starting to 
oscillate. 

Equations ()55II57|1 describe the curvaton model if we set Cm = Cr = jUl in which 
case ip represents the curvaton during its oscillating phase, with (^^in = Ca,in ^ Cm- By 
considering a massive curvaton we have |^ 

-H^ = -^ = -— (75) 
Pa 3 Pa 3 cr ■ 

For the last equality we have assumed that the field remains overdamped until the Hubble 
parameter falls below the curvaton mass, which is the case if Eq. (|7ij) is satisfied. Thus, 
from Eq. (j7!2|) we obtain 

C = r'^— = -raa — . (76) 

If we want to compare the efficiency of the inhomogeneous reheating with that of the curvaton 
model we must compare the efficiency parameter cto^ = — | to in Eq. (jSHI), which varies 
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according to the dependence of the mass and decay rate on the hght field 0. However, in 
the simplest case where F oc M oc 0, has opposite sign to and it is nearly twice as 
larger. 

We end this section with a comment. The physical situations discussed in this section 
using Eqs. ()55II57|) assume that the mass of must remain smaller than the Hubble parameter 
during the whole process of inhomogeneous reheating. However, this relation may be violated 
during the ip domination if ~ F and starts to oscillate during this period. This leads to 
a mix situation of curvaton/inhomogeneous reheating scenario, where the perturbation of 
is converted into a curvature perturbation via both the curvaton mechanism (the curvaton 
being 0) and the inhomogeneous reheating. Since the sign of the efficiency parameters a, 
can be different for the curvaton and inhomogeneous reheating models this may lead to 
a compensation between them. An example of this situation is given in The full 

calculation of the resulting curvature perturbation in the inhomogeneous mass-domination 
can be done starting from Eqs. ()39H43p . 



V. OBSERVATIONAL CONSTRAINT OF THE MODEL 

Here we discuss the observational predictions of the inhomogeneous reheating models: 
isocurvature perturbations and non-Gaussianities. 



A. Isocurvature perturbations in the mass-domination mechanism 

If the inflaton or the ^/'-particles decay into species out of equilibrium which remain 
decoupled from the radiation, we expect isocurvature perturbations to be present into these 
species. These can be correlated with the adiabatic perturbation. Here we consider the case 
where the perturbations left over from inflation are of the same order of magnitude as the 
perturbations produced during the inhomogeneous mass-domination and decay. We define 
the parameter n to quantify the relevance of the perturbations left over from infiation, 

" = = (77) 

where x is the inflaton, V^-^ is the derivative of the inflaton potential, and H is the Hubble 
parameter, all evaluated at horizon crossing during inflation. When k is order unity or 
smaller, perturbations from inflation are important. For chaotic inflation we have 

/tchaotic = 2ra0— -. (78) 

Using an inflaton fleld which is ~ lOMpi and an efficiency ra^j, ~ 0{1) we obtain ^chaotic ~ 
O.2Mpi/0 so that perturbations from infiation are important if the VEV of is sufficiently 
large, ~ Mpi. 

The following analysis holds also for the curvaton model, although the VEV of the curva- 
ton during infiation should remain smaller than the Planck mass [see Eq. fj74|) ] and generally 
K S> 1. Thus, in the curvaton scenario it is not likely that the infiaton and curvaton gen- 
erated perturbations are of the same order. The VEV of the light field does not need 
to satisfy this constraint. Indeed, as discussed in [l^, if one wants to avoid that at high 
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temperature the non-zero density of the ^'-particles generates a large thermal mass for - 
which would make (f) too heavy and would spoil the simplicity of the mechanism - we must 
require ~ M-p\. We are hence motivated to consider the possibility of k being small, at 
least for chaotic inflation. In this case perturbations from inflation may not be negligible 
with respect to perturbations from the light field and a mix of the two may survive. Isocur- 
vature perturbations in the curvaton model are discussed in and in [s^, although 

these grou ps considered a more general set of possibilities than what considered here and 
in IstI. liil they performed a numerical analysis of the CMB data to constrain the curva- 
ton model. Here however we consider a different possibility i.e., that the cold dark matter 
(CDM) is a relic left over both from the infiaton and the if) decay. 

We write the x and quantum perturbations as Qy^ = {H/ V2k^)a^ and = 
{H/\/'2k^)d^, where the d^'s are independent normalized Gaussian random variables, obey- 
ing (aQ,(k)a/3(k')) = SapSik — k'). After reheating, if the relic product of the infiaton is 
decoupled from the product of ip, its uniform curvature perturbation is conserved and given 
by 

Cproduct of X ~ 'y ~ ^X^X' C^^) 

From Eq. (jH^ the uniform curvature perturbation of the product of the ip decay can be 
written as ^ 

Cproduct of ip = -r^Qx jT^'t' ~ ^ '*^0)Cx' (^O) 

^,x ^ 

where the first term in the right hand side comes from the initial perturbation and the second 
from the fluctuations of 0. Note that k, according to its definition (f77|) . is scale dependent. 
We can write it as ^(k) = Ko(fc/fco)^"'^^ where kq is scale free and is a reference pivot 
scale. The spectral index of k. An, can be expressed in terms of the difference between the 
spectral indexes of and Cx; = fiif, — n^, but here it is considered as a free parameter. 

To simplify the discussion we completely neglect the baryons and we concentrate on 
the CDM isocurvature mode, which is due to the difference between the uniform curvature 
perturbations of the CDM and radiation (e.g., photons and neutrinos). We start by assuming 
that both the infiaton and the '0-particles may decay in CDM particles which are out of 
equilibrium at the temperature at which they are produced. We define / as the fraction of 
CDM, evaluated just before nucleosynthesis, which is left over from the decay of ip. The rest 
of the CDM, 1 — / , is a relic of the infiaton decay. Both the infiaton and the '0-particles may 
decay into radiation. The fraction of radiation that decays from ip is proportional to the 
value of flfj, at the decay, which we assume to be n<^ max defined in Sec. IV. B. We have seen 
there that this is very close to r. If r -C 1, ■0 remains negligible before decaying and cannot 
be responsible for the radiation ji^. However, tp can generate part of the CDM leading to 
an uncorrelated isocurvature mode. This is constrained by data: we have Jhq < 0.28 at 
95% confidence level. These bounds come from the numerical analysis of the WMAP data 
made in ji^], which assumes -0.72 < An < 1.11 and ko = 0.05 Mpc"\ The amount of 
CDM produced by the decay of can be important only if the ip generated perturbation 
is negligible. More interestingly, when r ~ 1, -0 comes to dominate the universe and the 
totality of the radiation comes from its decay product. In this case the adiabatic and CDM 
isocurvature perturbations are correlated. The intermediate case, that the radiation left 
over at nucleosynthesis is created both by the infiaton and by the -0 is not discussed here. 
Indeed, the radiation thermalizes and it is difficult to express its final perturbation in terms 
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of primordial perturbations. 

We thus concentrate on the case r ~ 1. Making use of the notation of i40.3 (see also j4ll|). 
the adiabatic mode is described by the comoving curvature perturbation in the radiation- 
dominated era, 7?.rad = — C7; and the isocurvature mode by iSrad = 3(Cc — C7); where C7 and Cc 
are the uniform curvature perturbations of the radiation (i.e., photons and neutrinos) and of 
the CDM, respectively. Well in the radiation era they are both constant and can be written 
as 

"^rad = (k-CL^ - a^)Cx^ (81) 

^rad = 3(1 - f)Ka^Cx- (82) 
According to ji^ we define the dimensionless cross correlation as 

I'^ol 



. ('^rad'^rad) 
COS A ~ 



k=ko \ I + K, 



,2 



(83) 



and the entropy-to-adiabatic ratio as 

/iso = 



('^rad) \ 



(^?ad), 



= 3(1-/) COS A. (84) 

k=ko 

These depend on two parameters, / and kq. For large values of kq the adiabatic perturbation 
is dominated by the perturbation of the light field and the modes are totally correlated 
(cos A = 1) as found in j^^. More generally the correlation is positive but can be small if 
the inflaton perturbation becomes important. 

The entropy-to-adiabatic ratio is constrained by data and cannot be too large. We can 
lower it by decreasing the amount of relic CDM left over from inflation (i.e., by sending / to 
1) or by decreasing the amplitude of the light field perturbation, i.e., the amplitude of kq. 
Since a full analysis of the constraints imposed by the data on this model is well beyond the 
scope of this paper, we just use the 95% confidence level bounds on the isocurvature mode 
coefficient /iso as a function of cos A for correlated perturbations as given in (Fig. 1 in 
this reference) and we show the bounds on / as a function of kq in Fig.|3J Totally correlated 
perturbations are allowed only if the CDM is almost completely created hj (/ = 1). 
However, even for k ^ 1, ~ 10% of CDM created by the inflaton is allowed. The signature 
of a correlated CDM isocurvature mode in the data can be the signal that part of the CDM 
has been created before the decay of ^. These constraints apply equally well to the curvaton 
model in the situation discussed here. 



B. Non-Gaussianities 

Until now we have assumed that the density perturbation of the light field depends 
linearly on Q(f,/4>, which we take to be a Gaussian variable. However, when the perturbation 
is comparable to the average value - which is the case for small values of - the 
non-linear term {Q^/(j))'^ can be important and lead to a non-Gaussian contribution in 



the spectrum of curvature perturbations [16|, |2 
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FIG. 4: Bounds on /, the fraction of CDM created by ip, as a function of kq = ra<^Q<^/((/)(^^) at 
95% confidence level. The filled region is excluded. 



The level of non-Gaussianity is conventionally specified by the non-linear parameter /nl 
(4^ li^j]. We can write the total uniform curvature perturbation as [i^ 

C = C. - IfNLCl (86) 

where (g represents the Gaussian contribution to (. Using Eq. (jH3j) the prediction for the 
inhomogeneous reheating is 

which is the same as for the curvaton scenario once the replacement = —2/3 of Eq. (|7S|) 
is used. However, Eq. (jH7|) has to be taken with caution: In order to precisely estimate the 
non linear parameter /nl generated by these models one has to study and solve the second 
order perturbation equations as done in ji^ and find the second order correction to Eq. (jU^ . 

If we use Eq. (jH7|) we see that less efficiency in the mechanism of conversion of per- 
turbations means more non-Gaussianities in the spectrum. If detected non-Gaussianity 
could be the smoking-gun of models where perturbations are produced with an 'inefficient' 
mechanism of conversion. The WMAP experiment has now put a limit on /nl corre- 
sponding to —58 < /nl < 134 at the 95% level which already excludes models with 
—0.006 < a^j, < 0.014. Planck will put a more sever constraint, |/nl| ^ 5 corresponding 
to \a^\ > 1/6. 

Going back to the inhomogeneous reheating, in which the mechanism of conversion of the 
density perturbation of into a curvature perturbation is due to the fluctuations of the decay 
rate and F oc we have ra,j) < 1/6 and the non-Gaussianity can be large, —5 < /nl < 0. In 
particular, if the non-relativistic species ip completely dominates the universe before decaying 
(r = 1) /nl = —5 a value which is right in the ball-park of Planck observations. 
However, the inhomogeneous mass-domination can be much more efficient. If M oc F oc 0, 
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ra<^ < 3/2 and thus —5/9 < /nl < 0. If the massive species dominates completely the 
universe before decaying (r = 1) we have /nl = — |; a much smaller value than the one 
estimate in and not observable by future planned experiments. 



VI. CONCLUSION 



In this paper we have studied the evolution of perturbations during a phase dominated 
by massive particles whose mass and decay rate can fluctuate in space and time. If the 
fluctuations are set by the VEV of a light scalar field overdamped during inflation, the 
isocurvature perturbation in the scalar field can be converted into a curvature perturbation, 
and this can be the main mechanism of generation of large scale perturbations for structure 
formation. We have derived a set of perturbation equations that can be used in full generality 
for any kind of dependence of the mass and decay rate on the light field. Making use of these 



perturbation equations we have recovered the results of [12, llll, UM inhomogeneous 
reheating with varying decay rate. We have also discussed the condition for the massive 
particles to dominate the universe before decaying. This condition does not depend on their 
mass, but depends on the annihilation cross section and decay rate. Standard model massive 
particles cannot dominate the universe. Furthermore, we have shown that when the mass 
of the massive particles is allowed to vary, the mechanism of conversion can be nine times 
more efficient. The final total curvature perturbation is C = — (l/6)5r/r — {A/3)6M/M. 
This is our main result. Finally, we have compared this with the curvaton model discussing 
differences and similarities. 

On the observational side we have discussed two possible signatures of the mass- 
domination mechanism: correlated adiabatic and isocurvature perturbations and non- 
Gaussianities. If present, a cold dark matter isocurvature perturbation provides some 
important information on the mechanism of generation of the dark matter and on the 
vacuum expectation values of the inflaton and light field during inflation. There are 
non-Gaussianities generated by this mechanism, which are x^- order to precisely 
compute them one has to study the evolution of second order perturbations. In the limit 
where /nl is large, by simply using Eq. (jHTj) the non-linear parameter is /nl = — 5/(6ra(^). 
When both the mass and the decay rate of the massive particles fluctuate, due to the high 
efficiency the non-Gaussianities can be much smaller than what is possibly observable, 
/nl < -5. 
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